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Observation of Coherently-Enhanced Tunable Narrow-Band Terahertz Transition 
Radiation from a Relativistic Sub-Picosecond Electron Bunch Train 
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We experimentally demonstrate the production of narrow-band (Sf/f ~ 20% at / ~ 0.5 THz) 
THz transition radiation with tunable frequency over [0.37, 0.86] THz. The radiation is produced as 
a train of sub-picosecond relativistic electron bunches transits at the vacuum-aluminum interface of 
an aluminum converter screen. We also show a possible application of modulated beams to extend 
the dynamical range of a popular bunch length diagnostic technique based on the spectral analysis 
of coherent radiation. 

PACS numbers: 29.27.-a, 41.85.-p, 41.75.Fr 



Techniques available so far to produce THz radiation 
are either limited (optical techniques do not easily pro- 
vide control of the THz pulse properties [H) or expen- 
sive (free-electron lasers require the operation of a large 
electron accelerator 0). Motivated by these facts, sev- 
eral groups are exploring alternative ways of producing 
THz radiation either based on low-energy (keVs) elec- 
tron beams 0-0] or compact (MeVs) accelerators @- 
[Toj . In this Letter we demonstrate the use of a 14-MeV 
pre-modulated electron beam to generate coherently- 
enhanced transition radiation (TR). Transition radiation 
is produced whenever a charged particle crosses the inter- 
face of two media with different dielectric constants [ll[ . 
The TR produced by a single particle is spectrally broad- 
band. The method we present is general and can be 
straightforwardly applied to any radiative electromag- 
netic processes jl2| . 
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FIG. 1: (color online) Bunch form factor S(u) associated with 
a train of Nb identical Gaussian bunches separated by T = 6at 
(a) and T = Wat (b) where at is the bunches' rms duration. 

In general, the spectral angular fluence emitted by 
a bunch of N 3> 1 electrons from any electromagnetic 
process is related to the single-electron spectral fluence, 
>n 1 ; " 1 ~ ^1 r * r -■- N 2 \S(u)\ 2 } where 

the bunch form 
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uj = 2irf (/ is the frequency) and S(uj 
factor (BFF), is the intensity-normalized Fourier trans- 
form of the normalized charge distribution S(t) [Uj]. The 
former equation assumes the bunch can be approximated 
as a line charge distribution and is practically valid as 
long as the rms bunch duration a t and transverse size 
a± satisfy a±_ <C cat/l where 7 is the Lorentz factor 
and c the velocity of light. When the BFF approaches 



unity, 
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N 2 and the radiation is termed as 



"coherent radiation" . Considering a series of Nb iden- 
tical bunches with normalized distribution A(t) we have 

S(t) = N^ 1 Y,nLi A (t + nT ) (where T is the period be- 
tween the bunches) giving | ^(oj) | 2 = £|A(w)| 2 where A(ui) 
is the Fourier transform of A(t) and the intra-bunch co- 
herence factor £ = iV~ 2 sin 2 {uN b T /2)/[sin 2 {uT /2)\ de- 
scribes the enhancement of radiation emission at the res- 
onant frequencies u) n = 2nn/T. The resonant enhance- 
ment is characterized by a full-width half-max (FWHM) 
value 8f ~ 2f 1 [ v / 3/(N 2 - IJJ/tt (for N b > 1). The BFFs 
computed for bunch train with different parameters are 
displayed in Fig. [TJ 

The production of a train of bunches with temporal 
separation T < 1 ps suitable for producing THz radia- 
tion is challenging. Most of the techniques developed and 
tested so far have severe limitations 0, H3, E3, EH • In the 
experiment presented in this Letter, we implemented an 
alternative scheme based on the exchange of phase space 
coordinates between the transverse and longitudinal de- 
grees of freedom 0, [I?} ■ 

The experiment was performed at the Fermilab's AO 
Photoinjector [l8[ ; see Fig. [21 In brief, electron bunches 
are generated via photoemission from a cesium telluride 
photocathode located on the back plate of a 1+1/2 cell 
radio-frequency (rf) cavity operating operating on the 
TM010 7r-mode at 1.3 GHz (the "rf gun"). The rf gun is 
surrounded by three solenoidal lenses (LI, L2, and L3) 
that control the beam's transverse size and divergence. 
The beam is then accelerated in a 1.3-GHz superconduct- 
ing rf cavity (the booster cavity) to ~ 14 MeV. Down- 
stream of the booster cavity, the 500-pC bunch is inter- 
cepted by a multislit mask consisting of 48- /im wide slits 
with 1-mm spacing thereby producing a transversely- 
segmented beam with total charge of ~ 15 pC. The beam 
is transported, with a set of quadrupole magnets, to the 
phase space exchange (PEX) beamline which consists of 
a liquid-Nitrogen-cooled deflecting cavity operating on 
the TMno-like 7r-mode at 3.9 GHz, located between two 
dispersive sections; see Fig. [5] The simulated trans- 
verse and longitudinal phase spaces evolution along the 
accelerator beamline is shown in the sequence of plots 
(b-g). The simulations, performed with the particle- in- 
cell package IMPACT [19( , illustrate how the transversely- 
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FIG. 2: (color online) Diagram of the AO photoinjector with 
the distances in mm (a). The "L", "Q" and "D" respec- 
tively refer to solenoidal lenses, quadrupole and dipole mag- 
nets. The density plots are snapshots of the centered trans- 
verse (Ax,Ap x ) (b-d) and longitudinal (Az,Ap z ) (e-g) phase 
spaces upstream (b,e) and downstream (c,f) of X3, and at 
X24 (d,g). The red traces are the associated horizontal (up- 
per row) and longitudinal (bottom row) projections. 



at fi ~ 0.26 THz of the beamline's transmission func- 
tion [see Fig. [3] (b)] is mainly a consequence of the fi- 
nite size of the X24 aluminum mirror and vacuum win- 
dow and was simulated using a vector-diffraction model 
for propagating the TR wave- front 20]. At high fre- 
quencies, the measured transmission for the z-cut quartz 
vacuum window was used pli] giving a high-frequency 
cutoff at f u ~ 2.5 THz. During our experiments, the 
beam was focused on the TR converter with the help of 
quadrupole magnets Q3, Q4 and Q5 to crj_ ~ 400 /ton. We 
therefore expect the beam transverse size to introduce a 
high-frequency cutoff at f s ~ 7c/[27rer^] ~ 3.5 THz be- 
yond the cutoff imposed by the z-cut window. The finite 
bandwidth A/ = [0.26, 2.5] THz of the detection system 
does not affect the measurements reported below but ulti- 
mately prevents us from resolving bunch structures below 
- 20 fim (70 fs). 

For current-modulated beams, the autocorrelation 
function is multipeaked and the peak separation provides 
the bunch-to-bunch separation within the bunch train. 
A Fourier transform of the autocorrelation function di- 
rectly provides the spectral power density associated with 
the detected radiation. The spectrum can then be used 
to infer the resonant frequency and the associated reso- 
nance width. An example of measured autocorrelation 
and computed peak separation are shown in Fig. [4] The 
bunch spacing is T = 2.3 ps and the spectrum, computed 
using a fast-Fourier transform (FFT) algorithm, exhibits 
the corresponding fundamental /i ~ 0.43 THz and sec- 
ond harmonic f% ~ 0.86 THz resonant frequencies. 



segmented beam upstream of the PEX beamline [Fig. [5] 
(c)] is converted into a longitudinally-segmented beam 
downstream of the PEX beamline [Fig. [5] (g)]. Varying 
the settings of the quadrupole magnets upstream of the 
PEX beamline [Ql and Q2 in Fig. [2] (a)] enables the 
control of the final longitudinal phase space parameters 
and consequently the bunch train parameters. Down- 
stream of the PEX beamline, coherently-enhanced tran- 
sition radiation (CTR) is generated as the beam transits 
the vacuum- aluminum interface of an aluminum TR con- 
verter labeled as X24 in Fig. [3] The radiation pulse is 
extracted through a z-cut crystal-quartz vacuum window 
and transported to a Michelson interferometer where the 
signal is recorded with a Helium-cooled bolometer; see 
Fig. [3] (a). The low- frequency suppression and cutoff 
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FIG. 3: (color online) THz beamline and autocorrelator (a) 
and corresponding transmission function (b). PI, P2 are off- 
axis parabolic mirrors, Ml, M2 are planar mirrors, M3 a mov- 
able mirror, and BS is an Inconel™-coated beam splitter. 
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FIG. 4: (color online) Measured autocorrelation function (a), 
corresponding peak locations (b) and deduced spectrum for 
the detected radiation (c). The numbers on plot (a) refer to 
the peak number. 

A series of autocorrelation functions was recorded for 
different settings of the upstream quadrupole magnets 
Ql and Q2. The measurements taken for Q2=l A and 
different values of Q2, are displayed in Fig. [5] (a) and il- 
lustrate the capability of the method to produce bunch 
trains with variable spacing thereby providing control 
over the resonant frequency of the TR; see Fig. [5] (b). 
In our experiment we were able to scan the frequency 
over /i £ [0.37, 0.86] THz while maintaining a relative 
bandwidth below 25%. During this measurement, a train 
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FIG. 5: (color online) Measured autocorrelation functions 
(with offset for clarity) for different current settings (indicated 
on the plot) of the quadrupole magnet Ql (a) and correspond- 
ing fundamental frequency and FWHM bandwidth (b). 



of 6 bunches was produced as inferred by observing the 
beam's transverse density downstream of X3. The dis- 
crepancy between the theoretical model described earlier 
[red line in Fig. [5] (b) with Nb = 6] and a linear regres- 
sion of the experimental data (green line) is attributed 
to the non-uniformity in charge and final duration spread 
for the bunches across the train. Such non uniformities 
among the bunches result in a larger FWHM bandwidth. 

Finally, the analysis of the CTR spectrum provides 
information on the bunch's temporal distribution [22j . 
Widely-used CTR-based diagnostics such as these 
are limited by the finite bandwidth of the detection 
system. Our system cannot measure "long" bunch with 
rms duration larger than ~ 1/// ~ 3 ps; see Fig. [3] 
Modulating the longitudinal density with a frequency 
within the bandwidth of the detection system allows for 
the measurement of these long bunches 23]. A direct 
demonstration of this feature appears in Fig. [6] where the 
autocorrelation functions of a long bunch were measured 



with and without intercepting the beam with the X3 
slits. For the nominal bunch most of the BFF content is 
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FIG. 6: (color online) Autocorrelation functions obtained 
with (blue) and without (red) inserting X3 slits (a) and cor- 
responding bunch form factor (b). The green dashed line in 
(a) represents the envelop of the multipeaked autocorrelation 
function (obtained with X3 slits inserted). 



at frequencies below /; [shown as shaded area in Fig. [6] 
(b)] thereby resulting in a autocorrelation function with 
strongly suppressed signal for times \t\ > 3 ps; see Fig. [5] 
(a) . In the case of the modulated beam the modulation 
leads to an enhancement of the BFF at ~ 0.43 THz so 
that information on the nominal bunch appears within 
the detectable bandwidth resulting in a modulated 
autocorrelation with its envelope representative of the 
nominal bunch shape. 
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